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■ Abstract 

a, 

We consider quantization of the Baierlein-Sharp- Wheeler form of the gravita- 
tional action, in which the lapse function is determined from the Hamiltonian con- 
straint. This action has a square root form, analogous to the actions of the rela- 
tivistic particle and Nambu string. We argue that path-integral quantization of the 
gravitational action should be based on a path integrand expfx/iS'] rather than the 
familiar Feynman expression exp[zS], and that unitarity requires integration over 
manifolds of both Euclidean and Lorentzian signature. We discuss the relation of 
this path integral to our previous considerations regarding the problem of time, and 
extend our approach to include fermions. 



1 Introduction 



Square-root Lagrangians are a feature of many field theories which are invariant 
under a time-reparametrization. The action of a relativistic particle 

S P = -m J dr ^-g^drX^drX" (1) 

and the action of the Nambu string 

S N = T J dadr yJ(d T x) 2 (d a x) 2 - (d T x ■ d a x) 2 (2) 

are familiar examples. Somewhat less familiar is the Baierlein-Sharp- Wheeler (BSW) 
form of the gravitational action 

Sbsw = -^Jd 4 x Jy/g 3RGv™(d tgij - 2N (i;j) )(d t g nm - 2N (n , m) ) (3) 

which is obtained from the standard ADM action, as reviewed below, by solving the 
Hamiltonian constraint. It is well known that for a relativistic particle moving in an 
arbitrary curved background, and for gravity in general, the corresponding quantum 
theory lacks a well-defined probability measure and time-evolution parameter 0, 

i, i- 

In this article we will propose a path-integral formulation of these " square- root" 
theories which is something of a departure from the standard Feynman expression. 
For one thing, the integrand of our path-integral will involve an unconventional 
phase: 

expfv^S] rather than exp^S] (4) 

Secondly, we will regularize the integration measure so as to uncover what we believe 
to be the true time-evolution parameter of the quantum theory. Third, we will find 
it necessary to sum over path-segments of both real and imaginary proper-time, 
i.e. over time-like and space-like trajectories in the case of the relativistic particle; 
Lorentzian and Euclidean signature manifolds in the case of gravity. It will be 
shown that the combination of the regularization, the unconventional phase, and 
the inclusion of imaginary proper-time segments, leads to a unitary evolution of 
states which corresponds, via the Ehrenfest principle, to the standard classical 
dynamics. 

In two previous articles [|J], || we have advocated a transfer-matrix approach 
to quantizing time reparametrization-invariant theories. The present article essen- 
tially presents the "real-time" version of our former "Euclidean" approach. Our 
previous work did not include fermion fields, which involve certain complications 
in our formulation. In this paper, we will show how the fermionic fields are also 
incorporated into our approach. 



2 Minisuperspace Actions 



We begin by considering simple quantum-mechanical theories with a time-reparametrization 
invariance, i.e. the "minisuperspace" models of the form 



S = J dt{p a d t q a -NH[p,q]} 



H = ^-G ab (q)p aPb + mV(q) (5) 
2m 

where it is assumed that the supermetric G ab has Lorentzian signature (— +++... +). 
The " square- root" form of the action is obtained by solving for p a in terms of the 
time-derivatives of the {q a }, i.e. 

= n sh = NG^ 

opa m 

=► Pa = ^G ab d t q b (6) 
and then solving the Hamiltonian constraint for the lapse function 

= J- G «\ q ) PaPb + mV (q) 
2m 

m „ „ „„ b 



2N 2 



G ab d t q a d t q b + mV(q) 



= sj-^d^dtq* (7) 

Substituting (g) and ffi) into the minisuperspace action then gives the square-root 
form 

S = - m Jdt ^J-2VG ab d t q a d t q b (8) 

For V — |, this is simply the action for a relativistic particle of mass m, moving in 
a background manifold with metric G ab . 

In non-relativistic quantum mechanics, a path-integral is constructed out of 
elementary integrals which evolve the wavefunction by a small time-interval e, i.e. 

ip(x',t + e) = J d D x /i e exp[iS[(x', t + e); (x, t)]/h]ip(x, t) 

= U £ ^(x',t) (9) 

where S[(x', t + e); (x, t)] is the action of a classical trajectory between the points x 
at time t and x' at time t + e. The measure \i € is chosen so that tp{x,t + e) — > ip(x, t) 
as e — ► 0. With this rule, one finds 

U e = exv[-iHe/h] + 0(e 2 ) (10) 



where H (the Hamiltonian) is an e-independent Hermitian operator. Taking the 
e — > limit, the evolution operator for finite times 

U M = lim(f/ e ) A ^ 

e—fO 

= exp[-iHAt/h] (11) 

is a unitary operator. Straightforward imitation of this construction doesn't work in 
the case of the square-root theories, due to the time-reparametrization invariance. 
Because of this invariance, the action of a classical trajectory between an initial 
point q and an end point q' is independent of the time parameters t and t + e which 
label those configurations, i.e. 

S[(q',t + e);(q,t)} = S[q\q} (12) 

The resulting operator U e defined from (|S|) would therefore be e-independent, and 
also in general non- unitary. 

Let us see if it is possible to recover an evolution operator of the form (JIOD for 
the square-root actions, by making a slight change to the construction shown in eq. 
(||). The modification is to multiply the action S[q',q] by an e-dependent complex 
constant c e 

ip(q',T + e) = J d D q fx e exp[c e S(q',q)}ip(q,T) 

= UM,t) (13) 

which is to be chosen such that U e is a unitary operator (up to order e) of the form 

U e = exp[-iEe/h] + 0(e 2 ) (14) 

where JE is an e-independent operator, hermitian in the measure fi € . Begin, for 
simplicity, with a "minisuperspace" action having V = | and G a b = Vaf, i-e. the 
action of a relativistic particle in flat D-dimensional Minkowski space. Let 

x 'n = x m _ a x m ( 15 ) 

so that 

ip(x',T + e) — J d D x fM e exp[— c e m\J — r] a bAx a Ax b ]^(x, r) (16) 

Comparing this expression to the corresponding expression for a free non-relativistic 
particle 



ip(x , t + e) = J d D x /J, e exp 



S ii Ax t Ax J 
m- 



-ieh) 



1 



motivates us to try 

C = -4= (18) 
y—ien 

with the understanding that the "time" -step e now has units of action, and that 
the branch of the square-root is chosen so that the exponent in (|i~6|) has a negative 
real part. This choice does not quite complete the definition of U e , as there is still 
a question of the range of the integral over x. Should this integral range over all 
possible x, or should it be restricted so that the path-segment Ax a is timelike? To 
resolve this issue, we will compute separately the contributions from timelike and 
spacelike intervals. 

Following the usual steps leading from the path integral to the Schrodinger 
equation, expand ip(x, r) in a Taylor series around x' 



1p(x', T 



J d D x fi e exp[— m^J— r] flv Ax u Ax^/ V—ieh] 



x 



^0',t) + -^Ax u + I - ,.. A.r'-'A.r" - 



2 dx'^dx' 



U e iP(x',t) 



In order that 



lim JJ f 



the measure must be chosen to be 

[i~ x = J d D x exp[— m^J— rj^Ax^Ax^/ \f—ieh} 
Changing variables Ax — > x, we then have 



1 



J d D x [ix^x v exp[— m^J— rj tJLV x^x u / \/ — ieh] 

Denote x M = {t, x} and r 2 = x ■ x. Then, on grounds of relativistic covariance, 
1 



(19) 
(20) 

(21) 
(22) 



1 + 



d D x /j,— — exp[-m^—T] fil/ x^x l/ / \/—iefi] 



d^dy + ... 



1 



lB d 2 + ... 



(23) 



2{D - 1) I A 
where d 2 = r]^ u d^d u , and 

I a = / d D x exp[— m^J— rj^x^x" / y/—ieh] 

Ib = J d D x r 2 exp[— m^— rj^ u x^x u / V—ieh] (24) 

We first evalute I a ; the second integral Ib will follow easily. Starting with 

I A = a J dtdr r D ~ 2 exp [— m V t 2 — r 2 / V—ieh] (25) 



where 

2n (D-l)/2 

° = WT (26) 

divide the integral over t into two contributions, one from timelike and one from 
spacelike paths: 

Ia = cr J drr D ~ 2 Ij dtexp[— mVt 2 — r 2 /\/—ieh\ 

+ dtexp[-mVr 2 -t 2 /Vieh]j (27) 



where the branches of the square roots ^/—ieh and Vieh are taken with positive 
real parts, to ensure convergence of the integrals. Next 

I A = a J drr D ~ 2 iyJ dtexp[-mVt 2 - r 2 /y/-ieh] J dy2y5[y 2 - (t 2 - r 2 )} 
+ jT dtexp[-mVr 2 - t 2 /Vieh] J™ dy 2yS[y 2 - (r 2 - t 2 )]| 
= a f°° drr D - 2 f°° dy { — L= e -™y/V^ + V^ e -my/^ Q{r _ A 

Jo Jo [ + t* v r — y J 

POO 

= a drr D - 2 {F 1 (r) + F 2 (r)} (28) 
Jo 



where 

Fi (r) = f°° dv M - e - m y/^^' 

F 2 (r) = f dy V - e-^/^ (29) 

Jo v r ~~ y 

It is easy to see that asymptotically, as r — > oo, both .Fi(r) and F 2 (r) go like 
1/r. But that implies 

roo 

timelike paths contribution — cr I drr D ~ 2 F\(r) is divergent 

Jo 

roo 

spacelike paths contribution = a drr D ~ 2 F 2 (r) is divergent (30) 

Jo 

This means that if we were to restrict the paths to only timelike, or only spacelike 
paths, then I a (and Ib) would be hopelessly divergent, and the evolution operator 
U e would be ill-defined. The remarkable thing, which we now show, is that the sum 
of the two contributions is actually finite. 

Let us deform the contour of ^/-integration for the integral defining -F\(r) in eq. 
(p9|). As it stands, it runs along the real axis from to oo. Deform it to run along 
the imaginary axis from to —ir, and then parallel to the real axis from — ir to 
oo. There are no poles or branch cuts in the way, so the deformation is permissible. 
Then 



F 1 (r) = \l " dy+ 1' dy\ !> e -™v/V^ (31) 



\Jr 2 + y 



Change variables, y — > —iy 



Fi(r) = - / dy . / e -"»y/V55_ f . / e ~^/^ (32) 



jo \fr 2 — y 1 Jr a/?" 2 — y 2 

Its not hard to see that it is the first integral which causes the divergence of the 
r-integration. Adding together F\ and F 2 , observe that the first integral in (|32|) 
exactly cancels F 2 , leaving an expression which decays exponentially as r increases 

F, + F 2 = - dy y^ e -my/V^n (33) 
Jr v r y 

The contour of this integral runs parallel to the (positive) imaginary axis. Now 
rotate the contour by 90 degrees, so that it runs along the positive real axis. Again, 
there are no poles or branch cuts in the way, and the integral is convergent along 
any contour intermediate between the initial contour, and the 90 degree rotated 
contour. This gives 



Fi + F 2 = i I dy r71 — 71 e 



y_ p ~ m y I 'v 'ieh 

\Jy 2 — r 2 



'ieh 

Inserting this result into (p8|), one finds 

/ drr D - x K x 
Jo 



trK^-^L) (34) 



1(J 



mr 



-_ 2f Jieh. Dr ,D + U„,D-h 



1(72 ( ^T } r( ^- )r( ^- } (35) 



It is trivial to repeat all the above steps for the Ib integral, and the result is 

POO l YTi r f 

I B = ia drr D+1 KA "" 
Jo 



ieh 



, D , Vieh D+2 D + 1 . 



= za2-(^)- +2 r(-— + l)r(— — + 1) (36) 
m 2 2 

Finally, inserting (|35D and ( |36|) into the expression for the evolution operator, eq. 

(EST), we obtain 



U e = l+teh { -^^- V ^d^ u + 0(e 2 " 
2m 2 

= exp[-iMe/h] + 0{e 2 ) (37) 

where 

M=-^(D + l)rj^d v (38) 



7 



The operator M is e-independent, and clearly Hermitian for inner products 

< HQ\*P2 >= J d D x ii e Tp*(x,T)Qip 2 (x,T) (39) 
The evolution operator 

U At ee lim([/ e ) AT / e 

= exp[-iMAr/h} (40) 

is therefore unitary, as in the usual path-integral approach for theories without a 
time-reparametrization invariance. 

It should be emphasized that the unitarity of our proposed evolution operator 
depends both on the choice of complex constant c e = 1/ ^/—ieh, and also on summa- 
tion over both timelike and spacelike path segments. A glance at equations (|35| ) and 



(RH) shows that the crucial factor of ie in ( |3"7| ) could only be obtained if a 1/ \/—ie 
factor multiplies the action in (|13|). Furthermore, the fmiteness of the result de- 
pends on keeping contributions to the integrand from both timelike and spacelike 
path-segments; the integral over either contribution separately is divergent. 

It is easy to generalize from the relativistic particle action to any minisuperspace 
square- root action of the form (H). First define the modified supermetric 



g ab = 2VG ab (41) 



S[q',q} = -mj^dt ^-Gabd t q a d t q b 



For Aq a small 



= -m\J-g ab Aq a Aq b (42) 

The measure is 

H-\<J) = (VeJ) D lim J exp[-mJ-g ab Aq"Aq>/V^idi] (43) 



eh) 1 

Now introduce Riemann normal coordinates £ a around the point q' a , which trans- 
forms the modified supermetric into the Minkowski metric Q ab = t] ab at the point q' 
(f = 0). Then 



, r dAq a 



b 



' Ia (44) 
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Inserting this measure, and c e = l/^—ieh, into eq. (|l3l), we have 



— f d D Aq^/\g(q> + Aq)\exp[-m\/-g ab Aq»Aq b /V^idi} 
xip(q / + Aq,r) 

j- J d D £ (1 - \n ab Ce + ...) expf-m^-^ee/v 7 ^] 



1 + (z^T)^ 2 -^ + °(« 2 » 



. (£> + i) afe a 2 

+ 2m 2 77 de^e 



or??/ 



^(g'(O^) 



(45) 



where TZ is the curvature scalar formed from the metric Q ab . Transforming back 
from Riemann normal coordinates, we have 



ip(q,r + e) 



ieh 



(D 



d 



2m 2 ^J\G~\ <^Z a 
exp[-iMe/h] + 0(e 2 )] tp{q, r) 



Q\Q 



lib 



d_ 

dq b 



ip(q,r) 



= U € ip(q,r) (46) 

As in the relativistic particle case, exp[— iJEe/h] is a unitary operator, where 

.,(£> + !) 1 



iE = -/V 



2m 2 



Q \ dq a 



oq° bm z 



(47) 



is obviously Hermitian in the measure jj, e of eq. (f44|). Taking the e — > limit, the 
wavefunction -^(g, r) satisfies a Schrodinger equation 

ihd T ip(q,T)=Mip(q,T) (48) 

The Schrodinger evolution equations (|47|) and ([48"D for time reparametrization 
invariant theories have been obtained by us previously, in refs. @] and f|, from a 
transfer matrix approach. The transfer matrix involves integration over a purely 
real integrand, but in our case the cost was not simply a Wick rotation of the 
evolution parameter r, but also a rather unnatural rotation of signature of the 
modified supermetric Q ab from Lorentzian to Euclidean. This rotation then had to 
be undone in postulating the Schrodinger equation (f48|). We have now seen that this 
supermetric signature rotation can be avoided, and unitary Schrodinger evolution 
is derived directly. In refs. || and || the correspondence of this evolution to the 
usual classical dynamics was also discussed. In the interest of completeness we will 



o 



briefly review this correspondence here, and refer the reader to the cited references 
for further details. 

The classical quantity M c i corresponding to the operator JE is obtained by re- 
placing derivatives with c-number momenta 

&ci{q a ,Pa} = ~» Pa] (49) 

which gives 

M d = (D + l) ^ GabpaPb (50) 
mv 

The Poisson bracket evolution equation 

d T Q = {Q,m cl } , M d = -e (51) 

is easily checked to be equivalent, up to a time-reparametrization, to the standard 
brackets 

d t Q = {Q,NH £ } , H £ = (52) 

where 

H £ = ^-^=G ab PaPb + VSmV (53) 

The parameter 8 is classically irrelevant, in the sense that it drops out of the Euler- 
Lagrange equations; the square-root action corresponding to 7i £ is 

S £ = -VSm J dt ^-2VG ab d t q a d t q b (54) 

Since £ only appears as a parameter multiplying the action, the fact that it drops 
out of the Euler-Lagrange equations is obvious. The same can be said for the mass 
of a relativistic particle in free-fall, the tension of the Nambu string, or Newton's 
constant in pure gravity. None of these parameters appears in the equations of 
motion at the classical level. 

Because of the classical equivalence of the Poisson bracket equations ( |5TD and 
(|5^), it is clear that the Schrodinger evolution (pEBD will obey an appropriate Ehren- 
fest principle, with certain quantum corrections due to the measure. The general 
solution i/j(q,r) can be expanded in terms of stationary states 

^(q,r)=J2Mq)e t£T/n (55) 
s 

where 

m<f> £ = -£<f> £ (56) 



1 n 



This T-independent equation can be rewritten in the form 



h 2 (D + 1) V d y _ x [7^ G a b d_ + h 2 (D + 1) 



S Am 2 



g\ dg a 



dq b 



£ Qm 2 



-vn + v 



<Pe(q) = (57) 



which is a Wheeler-DeWitt equation with a particular choice of operator-ordering 
and a (classically irrelevant) parameter S, which can be absorbed into a redefinition 
of either m or h. 

In the standard Dirac canonical quantization of actions of the form ([5]), the 
physical states must satisy a Wheeler-DeWitt equation of the form 



-ti 2 

2m 



'G 



ab 



d 2 



QqaQqb 



" +mV 



0(g) 







(5f 



where the quotation marks indicate an operator-ordering ambiguity. However, mul- 
tiplying the action by an arbitrary constant \f£ converts this constraint to 



n £ <i>(q) 



-h 2 d 2 

n G ab - - , " + V£mV 



2my/S dq a dq b 







(59) 



Because S is irrelevant at the classical level, there is no overriding reason that it 
should be a fixed parameter at the quantum level. In essence, our approach enlarges 
the space of physical states as the Hilbert space spanned by all states satisfying 



H £ (j>s{q) = o 



(60) 



and it is this enlargement of the space of states which enables us to obtain non- 
stationary states ip(q,r). Moreover, our approach fixes the operator-ordering, as 
seen in eq. (0), at least for the minisuperspace theories. Further discussion of 
these points may be found in ref. p. 

Returning to eq. (0), the path integral for square- root actions is now defined 
as the limit 



^(g/,r + Ar) = J Dq(r <r<T + Ar)e coS ^ T ^(q ,T 

Ar/e-l 

= lim / TT d D q n fi € (q n ) exp 

n=0 

lim(tg A ^%/,T ) 



Ar/e-l 



V-ieh m=0 



(61) 



where At = Tf — tq. We see that the time parameter emerges from a regularization 
of the path-integral measure: at fixed e, a regularized path between the initial point 



1 1 



go and the final point qf consists of n p = Ar/e path segments, each segment being a 
classical trajectory between intermediate points q n and q n+ \. The evolution param- 
eter is therefore a measure of the number of independent configurations (points) n p 
in the path joining q to qf, multiplied by the regularization interval, i.e. 

At = n p e (62) 

This is a quantum-mechanical time variable with no direct connection to, e.g., the 
proper time lapse. Nor is it an "intrinsic" time variable; all dynamical degrees of 
freedom are treated on the same footing and none is singled out as an evolution 
variable. Our evolution parameter is here identified as proportional to the number 
of "quantum steps" taken by the system in evolving from the initial to the final 
configuration. In this formulation the Green's functions are transitive, and the 
evolution of states is unitary. 



3 The BSW Action 



As in the minisuperspace case, the square-root form of the full gravitational action 
is derived from the first-order ADM action by solving for the lapse function. The 
ADM action for pure gravity is 



S = (fx \p lJ dtgij -NH- NiH 1 } 



n l = 



Jk 



ij^kl 



Vg 3 R 



IP 1 "* 

\9ik9jl + 9il9jk — 9ij9kl) 



2^9 



(63) 



where is the metric of a 3-manifold and 3 R is the corresponding scalar curvature. 
The momentum is related to the time-derivative of the metric by 

OH 



d t 9ij = N 



2K 2 NG ijk ip kl + N Mj + N j;i 
? G^ kl (d t9kl -2N m ) 



2k 2 N 

Solving the Hamiltonian constraint for the lapse function 



(64) 





N 



' G^ k \d t g l3 - 2N {i , j) ){d t g kl - 2N m ) - \^fg 3 R 



Ak 2 N 2 
1 



K 



G^ kl (d tgij - 2N {i , j) ){d t9kl - 2N m ) 



1/2 



(65) 
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and replacing the momenta in (|63f) by the expression (|64l) , with lapse (|65|) , gives 
the Baierlein-Sharp- Wheeler action |J 

5 BSl y = ~ j d A x WG^idtgij - 2N {l , j) )(d t g nm - 2N (n . m) ) (66) 

Before quantizing, it is convenient to fix the coordinate system by choosing shift 
functions Ni = 0. Then the corresponding supermomentum contraints 5S/5Ni = 
Ji l = are to be imposed as operator constraints on the space of physical states. 

It is straightforward to extend the BSW action to include non-gravitational 
bosonic fields. To compress indices somewhat, we introduce the notation 

{a = 1 - 6} i < j} 

q a (x) «-> gij(x) 

j p v(x) (i = j) 



Pa{X) 



2p lJ (x) (i<j) 



G ab (x) <- G^ nm (x) (67) 

and the non-gravitational fields are represented by q a (x) with indices a > 6. It is 
convenient to rescale all non-gravitational fields by an appropriate power of k so 
that all fields, and all components of the supermetric, are dimensionless. The action 
is 

S = J d 4 x[ Pa d t q a — NTt - NiTC] 

H = K 2 G ab PaPb + ^U (68) 

where ^ 

yfgU = ^^/g 3 R + non-gravitational contributions (69) 

Kj 

Setting the shift functions to zero and repeating the above steps of solving for the 
lapse, gives again a square-root action 

S = -- / d A x ^-^UG ab d t q-d t q b (70) 

The next step is to construct the evolution operator U e for the BSW action in the 
path-integral approach, following the procedure of the last section. The evolution 
operator is defined by 

V[j{x), r + e] = J Dq{x) ^e^^ 9tf(x)] + J d 3 x (J^) WW 

4/^ 3 ^^l^) Ag0(a;)Ag6(2/) + -- 
= ^(g , ,r) + [T + r 1 + r 2 ]+O(e 2 ) 

= UMq',t) (71) 



where the T n represent terms with n functional derivatives of \& and one power of 
e, and 

AS = -- f d 3 x J-^UG ah Aq-Aq* 

Aq a = q a -q' a (72) 
In order to obtain U £ , we need to evaluate 

< Aq^x^Aq*^) > = J D(Aq) (fi) Aq a (x 1 )Aq b (x 2 ) 



x exp 

where Q a b is the modified supermetric 



— / d 3 x(^) ^/-(g ab ) Aq-Aq b /V^Ji 



(73) 



g ah = ^-UG ab (74) 

and ()o indicates that the quantity in parenthesis is evaluated at q = q'. 

Clearly, < AqAq > is a highly singular quantity, and is only well-defined in the 
context of a regularization procedure. In the absence of a non-perturbative regulator 
which preserves the exact diffeomorphism invariance, we work with a naive lattice 
regulator in which the continuous degrees of freedom labeled by x are replaced by a 
discrete set, labeled by n, associated with regions of volume v n . We have in mind, 
e.g., a Regge-style discretization of a continuous 4- manifold into a fixed number N p 
of simplices of varying volume. As N p — > oo, the choice of v n is of course required to 
be irrelevant in computing the evolution operator, as long as the regions' volumes 
v n — > in this limit. As we will see, this requirement is not satisfied trivially or 
automatically. We take the naive lattice-continuum correspondences to be 



Aq a (x) <-> Aq a (n) 
f JV P 

/ d 3 x^/g <-> J2 v n 

J n=l 



5 f 6 \ \g[n) d 



5q a (x) \5q a (n) J R v n dq a (n) 

Dq <-> l[d D q(n) (75) 



With such a discretization, we have 

< Aq a (n)Aq\m) > = J \{d D q{j) fiAq a {n)Aq\ 



The supermetric G a b for the discretized degrees of freedom {q a (n)} still has 
Lorentzian signature, and we can follow the steps of the last section in integrating 
over the q a at each n. The result is 

1 



< Aq a (n)Aq b (m) >= teh(D + l) K 'Q ab —5 mn (77) 

v 



and we find for the T2 term 

2 V ^<9g a (n)<9g 6 (n) 
. t (£> + !) 2 v- ^ 1 ^ ab 1 d 2 

The term in braces has a simple continuum limit and, if this term were weighted 
by a volume factor v n , then the continuum limit would be simply 

2 J 5q a 5q b 

There is, however, no such v n weighting factor in the sum, which means that the 
contribution of each term at each position n 

1 w,; 1 f^y * (so) 



T 2 = ieh K ^ ' V / d 3 x U~ l G ab 1 — 1 — b ^ (wrong) (79) 



y/gU \6q a (n)J R \5q b (n) 

is weighted equally, regardless of the cell or simplex volume v n . As a consequence, 
even in the v n — > limit, the final answer for the state evolution would seem to 
depend on the distribution of volumes {v n }. 

Such regularization dependence never arises in ordinary quantum field theory. 
There may be other regularization issues, such as renormalization and anomalies, 
but certainly one doesn't encounter this kind of dependence on the distribution of 
cell volumes in computing the naive continuum limit of the Hamiltonian operator. 
Since the problem doesn't arise in ordinary quantum field theory, why does it come 
up here? The reason, of course, can be traced back to the square-root form of the 
action, which gives a factor of rather than a factor of l/v n , in the correlator 
(|T7j). The additional power of l/v n is the source of the (apparent) difficulty. There is 
only one way out, if the evolution operator is not to depend on the {v n } distribution: 
we must impose a constraint on the physical states $ such that the term j[8(\ ) is 
independent of the discretized position label n, at least in the v n — ► limit. In that 
case, we have 

Tl = fcS £+W iW t ) ( s ) , (anyn) (81) 

1 K 



The other terms To and T% are operator-ordering contributions which, in the ab- 
sence of an exact diffeomorphism-invariant regulator, will not be considered further 
here. Now absorbing the factor \{D + l)N p into a redefinition of e, and taking the 
continuum limit, we arrive at 
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Expanding 



v%, r) = e^<M<?) 
e 



(82) 



eq. 



requires that for each stationary state 
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~ n % Gab iAir- h 
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which is simply the Wheeler-DeWitt equation (up to operator-ordering contribu- 
tions), with an effective value of Planck's constant rescaled by 



ft - h 

heff ~ Te 
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Moreover, the Wheeler-DeWitt equation is consistent with, and in fact implies 
(via the Moncrief-Teitelboim interconnection theorem @), the supermomentum 
constraints 

W&e = (86) 

which are needed to compensate the gauge choice iVj = 0. Thus, each stationary 
state $g satisfies the usual constraint algebra of general relativity, with a rescaled 
value of Planck's constant. The Hilbert space of all physical states is spanned by 
the stationary states, with all possible values of £ . Finally, multiplying both sides 
of (p2|) by N^U, where N is an arbitrary function, integrating over space, and 
applying the supermomentum constraint (p6|), we obtain the equation of motion 
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d s x NG ab (-fr 
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where 
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and 



N(x) = mp- 



N(x) 



Jd*x' ^gNU(q) 
with mp an arbitrary parameter of dimension of mass. 



(89) 



The evolution equation (p7|) was obtained in ref. f| by a transfer matrix ap- 
proach, and shown to correspond to the usual classical evolution via the Ehrenfest 
principle. Here we have instead used the path integral to obtain a unitary evolution 
operator, as in the real-time Feynman approach, and avoided the signature rotation 
of the supermetric which was required in deriving the transfer matrix. The main 
point of this section is that, in performing the "real-time" path-integral of the BSW 
action 

*[g/(x),r + Ar] 
= J Dq(x,r)e coS ^ T H[q (x),r ] 

Ar/e-l r 1 Ar/e-1 

= lim / [] D( ln( x )Ve[qn}exp —= 

e ^ oj n= |y-»< 

xty[q (x),T Q ] 
= hm(U e ) A ^(q f (x),T ) 

= e- iMA ^[q f (x),T ] (90) 

it is necessary, as in the minisuperspace case, to integrate over all possible paths, 
including those for which the lapse function 

1/2 



S[q m+1 (x),q m (x)] 



N(x) 



1 



:G ab d t q a d t q b 



(91) 



4^ 2 vW 

is imaginary. Real-valued lapse functions correspond to Lorentzian 4-manifolds, 
imaginary values correspond to Euclidean signature. If the paths are restricted to 
real N(x) only, then we find that due to the Lorentzian signature of the supermetric, 
the integrals in eq. ( [T6|) are singular despite the regularization. 

The conclusion is that in order to obtain a unitary evolution of states, we are 
required to sum over 4-manifolds of both Lorentzian and Euclidean signature, and in 
general over manifolds which may be Lorentzian in some regions, and Euclidean in 
others. This raises the obvious question of why spacetime seems to have Lorentzian 
signature, rather than Euclidean or mixed signature. 

The question "why is spacetime Lorentzian?" can be raised already at the level 
of classical general relativity. Einstein's equations themselves do not specify a 
choice of metric signature; there are Lorentzian solutions to these equations, and 
there are Riemannian solutions. Recently, solutions to the Einstein equations in 
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which part of the manifold is Riemannian (Euclidean signature) and the rest is 
Lorentzian have been studied IIIBl; it is conceivable that solutions of this kind are 



relevant to the very early Universe. In any case, the signature of a manifold solving 
the Einstein equations is determined in general from initial conditions {gij,p^} 
satisfying the appropriate constraints. A given initial 3-manifold may trace out 
either a Lorentzian or Riemannian 4-manifold, depending on the initial choice of 
conjugate momenta. 

The dependence of lapse on initial conditions applies also to the quantum theory. 
The general solution of the evolution equation ( fI8|) for the "relativistic particle" 
example, with M given in eq. (|38D , is 



d 4 p f{p) exp 



i f( D + 1 ) 2 



(92) 



It is easy to see that 



< x M >=< x^ L > 4 



(D+l) 
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< jr > r 



(93) 



(recall that r has units of action). So long as f{p) = for p 2 > 0, the expectation 
value of position follows a timelike path. We would expect the same situation in 
quantum gravity, for the same reason, namely, the Ehrenfest principle. If the initial 
"wavefunction of the Universe" ty[q a (x), To] has expectation values which are peaked 
around some (equivalence class of) configurations and momenta {q,p}o, then the 
wavefunction tends to remain peaked in the neighborhood of a classical manifold 
which solves the Einstein equations for this initial data. Thus, despite the fact that 
the path integral sums over Lorentzian and Euclidean manifolds, the probability 
density can still be sharply peaked at one or the other signature. 

Obviously these remarks do not answer the question "why is spacetime Lorentzian?" . 
but only replace it with another question about initial conditions. For an attempt 
to explain the preference for Lorentzian signature (in the context of non-time- 
parametrized theories) from an analysis of an effective "signature potential", see 
ref. H. 



4 Quantum Theory in Curved Spacetime 

Associated with the problem of time in quantum gravity is a "problem of state." 
Let us return, for a moment, to the standard formulation of canonical quantum 
gravity, which in our language is a restriction to a single value of £, and let H be 
the Wheeler-DeWitt Hamiltonian. Suppose a physical state ^ is an eigenstate of 
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an observable Q; this means that Qfy must also be a physical state. But then 



H(QV) = [H,Q]* = 



(94) 



which is not true, in general, unless [H, Q] vanishes weakly. It is then problematic 
to construct physical states which are approximate eigenstates of, e.g., 3-geometry, 
or the position of the hands of a clock. 

In this section we show how to construct physical states which are sharply 
peaked around a given 3-geometry and extrinsic curvature. Treating the metric 
degrees of freedom semiclassically, the dynamics of the other degrees of freedom 
approximates the standard quantum theory on a curved background. Of course, 
the WKB treatment can be extended to any other degrees of freedom (such as the 
hands of a clock) which behave more or less classically. 

We recall that our path integral leads, in the end, to the following solution for 
the evolution of physical states: 



y[q,T] = ^c[£,a]<S>£ tCt [q]e 



iSr/h 



(95) 



where 

E<5> s ,a = -Z®s, a (96) 

and where the subscript a is meant to distinguish between different solutions of 
(j96|). As discussed above, eq. (|96|) is a one-parameter (£) class of Wheeler-DeWitt 
equations 



f? 5 2 
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(97) 



each of which can be treated by WKB methods. To get the quantum-theory-in- 
curved-spacetime limit, we follow the approach of Banks ||, treating the metric 
semiclassically, and expanding in powers of k 2 (back-reaction of matter on metric 
will be ignored; it can presumably be dealt with following the approach of ref. |10| ). 
Thus, write fl9"6T) in the form 



£ 



- k Gijkr 



t 2 1 ™ 







where l-i m is the Hamiltonian density for the non-gravitational fields, denoted 
We then make the WKB ansatz 



(9£ 



exp [iV£S[g,g]/K 2 h] p s vv [gW m [4>] 



(99) 



where S[g,g] is Hamilton's principal function (the action of a 4-manifold solving 
the Einstein equations, bounded by the three manifolds with metric g^ and gr^-); it 
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satisfies the Einstein-Hamilton- Jacobi equation in both arguments 

5S 5S 

ogij og k i 

G ^w^~^ m = (100) 

The functional Pyy[g] * s the Van-Vleck determinant, while i/'m is a solution of the 
Tomonaga-Schwinger equation for quantum theory on a curved spacetime back- 
ground 

foeff * = Ktf* (101) 



where 

Kff = -j= (102) 

is the effective value of Planck's constant, and T(x;g,g) is a functional of the 
background spacetime defined by 

61 dgij dg kl 

Up to this point, we have simply repeated the analysis of ref. || 

However, the semiclassical approach to recovering ordinary quantum field theory, 
as outlined above, is subject to the following objection: Although the part of the 
wavefunction involving the non-gravitational fields obeys a Tomonaga-Schwinger 
equation, the metric gij, on which the " many- fingered" time parameter T(x;g,g) 
depends, is still a dynamical degree of freedom, and there is no physical state 
satisfying the Wheeler-DeWitt equation which has a probability distri- 
bution peaked at a particular 3-geometry g^, i.e the wavefunction is not 
peaked on any particular time-slice of a 4- manifold. In fact, the squared- modulus 
of the leading term in the WKB approach, i.e. 

exp [iV£S[g,g]/K 2 h] | 2 = 1 (104) 

has no dependence on g^ at all. The best one can do in the standard formulation 
(that is, using only a single value of £ ) is to superimpose WKB solutions 

®eM = J Dg^Flg^ea 

Dg l3 f[g l3 ] exp [i{V£S[g, g\l k 2 + 9[g]}/h}} p £ vv ^ £ m (105) 



on 



where f[g] is a real functional peaked (modulo diffeomorphisms) at a particular 
3-geometry g 0i j, and we define 



Po 




(106) 
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where 9[g] is the phase of the smearing functional F[g]. As shown many years 
ago by Gerlach pA| , this superposition is still not peaked at any one 3-geometry, 
but rather on all three-geometries which are spacelike slices of a certain 4-manifold, 
satisfying Einstein's equations with initial data {goij,Po }• Thus there is no physical 
state, and no subspace of physical states, which would correspond to an eigenstate 
of a non-stationary observable (such as the three-geometry, or the fields on a given 
three-geometry) . 

It is at this point that we make use of the freedom, inherent in our formulation, 
to superimpose states of different £, and write 



d£Dg F[g, £} exp 



{£T + VSS[g } g}/K 2 } 



^[<j>,T{x;g,g )} J d£DgF[g,£]ex P ^-{£r + V£ S[g,g]/K 2 } 
^[0,T(x^,^ o )]4 9o ' po) [^,r] 
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(107) 



where it is assumed that F\g, £] is sharply peaked around £ = £ , g^ = g 0i j, and ptf 
is defined as in eq. (|106|) above. The ip m factor can be pulled outside the integral, 
on the grounds that its variation with £ and g is much less than that of the smearing 
function F[g,£], and the exp[i\/£S/K 2 h] factor. Now consider the leading WKB 
term 



4 5 °' P0 W^] = / dSDy^ exp y-{£ T + V£S[g,g]/K 2 + m} 
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This wavefunction will be peaked at configurations where the phase in the in- 
tegrand is stationary, with respect to small variations in 7^- and £ around goij and 
£ , respectively. In other words, the wavefunction is peaked at metrics gr^-, at time 
r, such that 
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V£o ( SS[g,g]\ 



(109) 



\g=go 



The second of these two equations is satisfied by the metric gij of any time-slice 
of a 4-manifold, satisfying Einstein's equations with initial data {goij,Po }• The 
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first equation requires that the action of the 4-manifold between the initial slice 
goij and the given slice is proportional to the time-parameter r. Now consider a 
foliation of the given 4-manifold parametrized by some variable xq, with g^ = g^ij 
at Xq = 0. Hamilton's principal function S[g,go] is monotonic in xq, which means 
that S[g,g ] = only for g^ = g 0i j. It follows that, at r = 0, eq. ( |109| ) gives us 

= S[g, go] =>- gij = goij (modulo diffeomorphisms) (HO) 

As a consequence, at r = 0, the wavefunction §p°' Po \g, r = 0] is peaked at g^ = g^j 
(modulo diffeomorphisms). Thus, from the definition of the many-fingered time 
variable, where T(x;go,g ) = 0, 

*foy>,T = 0] = ^°' Po) [g,r = 0] x ^,T(x;g,g )} 

« 4^)[^r = O]x^[0,T = O] (111) 



The importance of eq. (|1 1 1|) is that there exists, in our formulation, a class of 



states where the metric (and extrinsic curvature) is sharply peaked around a given 
geometry g 0i j (and Pq), and where the state factorizes into a wavefunction ($f) 
suppressing fluctuations away from the given 3-geometry, and a wavefunction (ip m ) 
describing the state of the non-gravitational fields on that 3-geometry. Such states 
can be fairly described as eigenstates of non-stationary observables; these eigen- 
states are impossible to construct, as physical states, in the standard formulation 
of canonical quantum gravity. 

Finally, we consider transition probabilities. Take an initial state of the form 

^ n [^,0] = $^^^,O]x^[0,O] (112) 

and a final state of similar form 

Vfha, 0] = ^?' P ' o) [g, 0] x 0] (113) 

where the smearing function F' is peaked around some time slice {g' ,p' ) of the 
classical 4-geometry specified by the initial data (<7o>Po)- The transition probability 
for — ► tyf after a time r is given by the factorized expression 

Pin-,f(r) = \<V f \e- iMT / h \% n >\ 2 

= | < * f \* in (T) > | 2 

= \<^ p ' \gM^ ' P0) [9,r]>\ 2 

x |<^t o [0,O]iC°[0,T(x;^^o)]>| 2 (H4) 
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The first of these factors 



<^' p '%,0)\^°' po) [9,r}> 



(115) 



gives the probability, after a time r, to be on the time-slice described by (g' ,p' ) (up 
to a certain uncertainty, specified by the smearing function F'). The second factor 



is the quantum- field-theory-in-curved-spacetime result; it gives the probability for 
a transition from an initial state ip m of quantum fields on the time-slice go, to the 
state i/j' m on the later time-slice g' . Both the initial and final 3-manifolds are time- 
slices of the same 4-manifold specified by the initial data {go, po}, and the state 
il> m [<t>, T) evolves according to the Tomonaga-Schwinger equation (|101|) . 

In this way, we see how approximate eigenstates of geometry and extrinsic curva- 
ture may be constructed, and how the standard formalism of quantum field theory 
in curved spacetime emerges. We will not attempt to go further and discuss the 
problem of measurement in this context, apart from noting that any of the standard 
"realistic" approaches that have been applied to non-relativistic quantum mechan- 
ics, e.g. many-universes, decoherence, or Bohm's theory, can be applied in our 
formulation as well. 

5 Inclusion of Fermions 

We have so far assumed that the canonical momenta p a appear quadratically in the 
Hamiltonian, with indices contracted by the supermetric. The Hamiltonian of a set 
of Dirac fields, on the other hand, is linear in the fermionic momenta, and it is not 
immediately obvious how such fields are incorporated into our approach. 

In our previous work || we found two independent methods for determining 
the M operator. The "undetermined constant" method was based on the trivial 
observation that the actions S and 5" = const, x S are equivalent at the classical 
level; this leads to the fact that the ratio of the kinetic and potential terms of the 
Hamiltonian (which is the M functional), is indeterminate at the classical level. The 
second method, leading to the same quantum theory, is the transfer matrix method, 
whose "real-time" or Feynman version was presented in the preceding sections. We 
will now apply both methods to obtain the M operator for gravity coupled to a 
Dirac field. 

The action for the Einstein-Dirac system is expressed in terms of the fermion 
field %fj and tetrad field e a ^ as 




(116) 




(117) 



where _D M is the usual covariant spinor derivative. The extension of the canonical 
ADM formalism to this system was worked out in ref. ||12||, for the "time-gauge" 



e ° i = (* = 1,2,3) (118) 

In this gauge, the Einstein-Dirac action expressed in terms of canonical momenta 
has the form 

S = J d 4 x [ Pa q a + Tfiptp - (NH + NiW + e ij Mij)] (119) 

where the q a are the triad fields e£(x), and 

H = K 2 G ab p a p b + ^U + H^ 

= n^KtP = 7V7 [el7 a A ~ (120) 

The first-class constraints are 

U = W = Mij = (121) 

where the supermomenta 7i l and the generators of local frame rotations Mij are 
linear in the momenta. In addition there are 2nd-class constraints, some of which 
are associated with the time-gauge condition (|118| ), and also which relate to ?/>:[] 



7T„ 



iy/gifa° (122) 



The 2nd-class constraints are handled, according to the Dirac procedure, by replac- 
ing Poisson brackets by Dirac brackets. The explicit form of all constraints in terms 
of the canonical variables, and other details, may be found in ref. |T2] ]. 

Now consider an alternative action S' ED which differs from Sed only by a mul- 
tiplicative constant 

S' ED = VSS ED (123) 

Obviously, the equations of motion derived from S' ED are identical to those derived 
from Sed- The constant £ is therefore irrelevant at the classical level. In going to 
the canonical formulation, however, we find that 

S' = J d 4 x [ Pa q a + vr^ - (NH £ + NiW + e lj M t 

n £ = -^=K 2 G ab Pa p b + H^ + Vs^u 

with 2nd class constraints enforcing 

tt^ = iVSy/g^ (124) 



1 We take right derivatives with respect to ip. 
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Define 



d 3 x (NH £ + NiT-C + e ij Mi, 



(125) 



and consider a field configuration {ef(x, t),ip(x, t), ip(x, t)} which solves the Hamil- 
tonian equations of motion derived from H , for some given value of £. Then it is 
clear that this configuration is a solution for any other value of £ , since the classical 
orbits in configuration space (i.e. solutions of the Euler-Lagrange equations) are 
independent of £. In general then, the Dirac bracket equation of motion Q 



d t F = {F,H £ } D 



(126) 



supplemented by the first class constraints 



H £ 



for any N, N h e* 



n l = Mi. 
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generates a set of orbits in configuration space which is independent of £ . In this 
sense £ is "classically irrelevant." 

Now observe that the constraint H = can be written 



E 



-£ for any N, N h e* 
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where E is defined implicitly by 
E = f d 3 x 



, j N(K 2 G ab PaPb + V^EH^ 

J d 3 x' ^fgNU in,, 



+ —(N l H i + e ij M ij ) 



(129) 



and where m p is an arbitrary parameter with dimensions of mass. From this defi- 
nition, it is straightforward to show that, for any functional F = F[q,p, if), tt^\, the 
Poisson bracket with E is related to the corresponding Poisson bracket with H via 



m p {F,E} 



d 3 x 



1 + 



/ d 3 Xl NH^ 



2V£ Jd 3 x 2 N^/gU. 
+N l {F,H t } + e^{F } M i3 } 

{F,H £ } N ^ 



\f~£m p N 
Jd 3 x 3 N^gU 



{F,H £ } 



where 



N(N) = 



1 + 



fd 3 x 1 NH, 



2y/£Jd 3 x 2 N^/gU 



\f£m p N 
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2 Of course, {F, H £ }d ss {F,H £ }, since Tt = is a first-class constraint. 
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Eq. ( p.30| ) is derived by simply carrying out the functional derivatives contained 
in the Poisson brackets shown, and applying the constraint ( [L28|) . Then, since the 



Dirac bracket {F, M}d is linear in Poisson brackets {..,JE}, eq. (|130|) implies 

m p {F,m} D = {F,H e }%-*X (132) 

Defining r = m p t, this demonstrates the equivalence of 

d t O = {0,H £ } D & d r O = {0,M} D (133) 

up to a time-reparametrization, expressed by N — > N. Note that N(x) and const, x 
N(x) have the same N. 

We now quantize by replacing Dirac brackets with commutators (in the case 
of bosonic fields), and anticommutators (in the case of fermionic fields). Time 
evolution of states is given by the Schrodinger equation 

ihd T ^[q, $,T]=Eti![q,ij), r] (134) 

with the general solution 

*[q,#,T] = Y,ae*e[q,'$\e ieT/n (135) 

£ 

where 

4> = g^ip (136) 
and satisfies a Wheeler-DeWitt equation 



£ 5q a 5q b </£5^ V 



$ £ = (137) 



where the operator-ordering remains to be specified. Note that, as in the purely 
bosonic case, the classically irrelevant constant £ can be absorbed into a redefinition 

of n 



h eff = — (138) 



This concludes the first, "undetermined constant" method for finding the M oper- 
ator. 

Next we apply the path-integral approach, following as closely as possible the 
procedure of the previous section for the purely bosonic case. Since the gener- 
alized BSW action for gravity + fermions will contain a factor of Ti,^ inside the 
square- root, our strategy will be to expand the path-integrand in powers of 
and evaluate the relevant expressions to some finite order (first order, in this arti- 
cle). These expressions can then be compared, order by order in 7i^, with results 
of the "undetermined constant" method above. 



We again set JVj = 0, and also e y = 0, which is to be compensated by imposing 
the corresponding physical state constraints 



m = MijV = 



Solving for the bosonic momenta in terms of the time-derivatives 
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2k 2 N 

inserting into the Hamiltonian constraint 
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G ab d t q b 
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and solving for the lapse function 

1 



G ab d t q a d t q b + y/gU + H^ = 
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N = — 



nl/2 



-G ab dq a dq b 



(142) 



2k [ (y/gU + H*) 
we arrive at a square-root action 

S = Jd 4 x {ijfh ^ ~ l^-(VgU + H^G ab d t q"d t q b } (143) 
The corresponding path-integral is 

for', t + At] = J DqD^Di, fx(q, ^, ?) exp c 1 d 4 x {iy/ff^d^ 

-- ^-(VoU + H^G ab d t q a d t q b }} *[qJ,T] (144) 

where 



H^ = i^fK^ 



(145) 



and Co represents the e — > continuum limit of the regularization-dependent con- 
stant c e oc \/yJ—ieh. Now expand the exponential to first order in H$ 

V[q','4>',T + At] 
= j DqD^Dil) 

-±J d 3 xdr -J^-^-^UG ab d T q-d T q b 

V[q,4>,T 



X 



x exp 

where S is the bosonic action 



c o ( So + / d 3 xdr iy/g^^fd^ 



S = -- / d 3 xdr^-^UG ab d T q a d T q b 



(146) 



(147) 
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We now regularize the path- integral according to the lattice prescription (f75). 
The bosonic part of this path-integral, based on the action S , leads to the operator 
evolution 



(148) 



where 



d 2 



dq a (n)dq b (n) 



and we have defined 
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dq a (n) 



(149) 
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Then to zeroth-order in Tl^, we can identify the derivatives d T q in the term pro- 
portional to Tirj} as proportional to the bosonic momentum operators, according 
to 
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dPa(n) 
(D + 1)k 2 



Uv r , 



-G ab p b (n) 



where we have introduced 



Pain) = —Pain) 



so that, as operators, using eq. ( |150| ), 
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The regularized path-integral, to first-order in TC^ is now 
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V[q' + Aq,4>' + Ai>,r] 
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where 
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Carrying out the integrals over Aq, Aijj, ijj, we find 

1 1 



r+e 



n 2 G ab PaPh 



-K 2 G ab PaPb ( JT 



Kip 
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(156) 



where the square-root operator is defined via spectral analysis. 

Once again, the requirement that the state evolution is independent, in the 
continuum limit, of the choice of {v n }, implies that the term in braces is the same 
in each cell (simplex) n. Therefore, in the continuum limit, 



K 2 G ab PaPb + 



-n 2 G ab p a p b 5 ~ 
, in — -Kw 



(all x) (157) 



where the divergent factor N p (D + l)/2 has been absorbed into a rescaling of r. 
For stationary states we have 



n 2 G ab p a p b + 



-K 2 G ab p a p b 5 ~ 
— th-jK-ip + S^/gU 



= 



(158) 



Then, since 



-K 2 G ab p a p b 



= VS&s + order 7-L, corrections 



(159) 



it follows that, up to first order in Ti^, we recover the same stationary state equation 

h 5 



h 2 fi 2 
£ 5q a 5q b 



+i 7tM K ^ + V ~ 9U . 



=o 



(160) 



that was obtained (eq. ( p. 37] )) from the "undetermined constant" approach. As in 
the purely bosonic case, the imposition of TC £ &£ = for all £ at every point x 
is consistent with, and in fact implies (c.f. ref. [0]), the other required first-class 
constraints 7i 1 ^! = = 0, up to the usual operator-ordering issues. 



6 Conclusions 

For ordinary quantum theories without time-parametrization, the regularized path 
integral is expressed as a product of integrals, each of which evolves the state 
function unitarily over a very small time interval. In this article we have examined 
whether such a construction can be applied to theories with square-root, time- 
reparametrization invariant actions. Our result is that unitarity requires: i) an 
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unconventional phase in the path-integrand; and ii) summation over configurations 
of both real and imaginary proper-time lapse. In the case of quantum gravity, 
the second requirement means that path-integration must run over manifolds of 
Lorentzian, Euclidean and, in general, mixed signature. We have also shown how 
the formalism extends to fermionic actions. 

Unitarity, of course, refers to evolution in a certain time parameter. In our 
formulation, the time parameter is simply a measure of the number of integrations 
in the (regulated) path-integral, evolving an initial state to a later state. This 
"quantum time" parameter is neither a geometrical quantity (such as a proper time 
lapse), nor a dynamical variable (such as the extrinsic curvature). It is, instead, 
a parameter which is intrinsic to the path-integral measure. The connection to 
classical dynamics is established via an Ehrenfest principle. 

In the standard canonical formulation of quantum gravity, the physical states 
are solutions of a Wheeler DeWitt equation Ti^ = 0. In contrast, an outcome of our 
formulation is that physical states belong to a Hilbert space which is spanned by 
the solutions of a family of Wheeler-DeWitt equations Ti <$>g = 0, which are distin- 
guished by having different effective values of Planck's constant h e ff = h/y/S. As 
discussed in section 4, a superposition of states with varying S (or ft e //) allows us to 
construct, at the semiclassical level, physical states whose amplitudes are peaked at 
particular 3-geometries and extrinsic curvatures. The width of the peak, in super- 
space, is inversely proportional to the dispersion AS. Projection operators formed 
from such states, and linear combinations of those projection operators, belong to 
the physical observables of the theory. It is worth noting that the stationary states 
(i.e. solutions of a Wheeler-DeWitt equation with a fixed value of S) can never be 
peaked around any one 3-geometry. At best, in the WKB limit, a stationary state 
is peaked at every possible spacelike slice of some 4-manifold satisfying the Einstein 
equations. 

If our view is correct, then the phenomenological value of Planck's constant is 
the mean value of a dynamical quantity, having a finite uncertainty of quantum 
origin. How large this uncertainty might be, and whether there could conceivably 
be testable consequences, are interesting issues for further investigation. 

Acknowledgements 

We would like to thank Maurizio Martellini for some very stimulating discus- 
sions. J.G. is also happy to acknowledge the hospitality of the Lawrence Berkeley 
Laboratory, and the Niels Bohr Institute. J.G.'s research is supported in part by 



the U.S. Dept. of Energy, under Grant No. DE-FG03-92ER40711. A.C.'s research 
is supported by an EEC fellowship in the 'Human Capital and Mobility' program, 
under contract No. ERBCHBICT930313. 



si 



References 



[1] K. Kuchar, in Quantum Gravity 2: A second Oxford symposium, ed. C. Isham, 
R. Penrose, and D. Sciama (Oxford University Press, Oxford, 1981). 

[2] K. Kuchar, Time and interpretations of quantum gravity, Proceedings of the 
4th Canadian Conference on General Relativity and Astrophysics, ed. G. Kun- 
statter et. al. (World Scientific, Singapore, 1992). 

[3] C. J. Isham, Imperial College preprint IMPERIAL-TP-9 1-92-25, 
archive: |gr-qc/9210011| . 

[4] J. Greensite, Phys. Rev. D49 (1994) 930. 

[5] A. Carlini and J. Greensite, Niels Bohr Institute preprint NBI-HE-94-31, 
archive: |gr-qc/9406044 . 

[6] R.F. Baierlein, D.H. Sharp, and J.A. Wheeler, Phys. Rev. 126 (1962) 1864; 
J. A. Wheeler, in Relativity, Groups, and Topology, ed. C. DeWitt and B. 
DeWitt (Gordon and Breach, New York, 1964). 

[7] V. Moncrief and C. Teitelboim, Phys. Rev. D6 (1972) 966. 

[8] A. Carlini and J. Greensite, Phys. Rev. D49 (1993) 866; 
J. Greensite, Phys. Lett. B300 (1993) 34. 

[9] T. Banks, Nucl. Phys. B249 (1985) 332. 

[10] R. Brout and G. Venturi, Phys. Rev. D39 (1989) 2436. 

[11] U. Gerlach, Phys. Rev. 177 (1969) 1929. 

[12] J.E. Nelson and C. Teitelboim, Ann. Phys. 116 (1978) 86. 

[13] G.F.R. Ellis, A. Sumeruk, D. Coule and C. Hellaby, Class. Quantum Grav. 9 
(1992) 1535; 

T. Dereli and R. Tucker, Class. Quantum Grav. 10 (1993) 365; 
T. Dray, C. Manogue and R. Tucker, Phys. Rev. D48, (1993) 2587; 
S.A. Hayward, Class. Quantum Grav. 9 (1992) 1851. 



10 



